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ABSTRACT 

Very-long baseline interferometric observations have resolved structure on scales of only a few 
Schwarzschild radii around the supermassive black holes at the centers of our Galaxy and M87. In the 
near future, such observations are expected to image the shadows of these black holes together with 
a bright and narrow ring surrounding their shadows. For a Kerr black hole, the shape of this photon 
ring is nearly circular unless the black hole spins very rapidly. Whether or not, however, astrophysical 
black holes are truly described by the Kerr metric as encapsulated in the no-hair theorem still remains 
an untested assumption. For black holes that differ from Kerr black holes, photon rings have been 
shown numerically to be asymmetric for small to intermediate spins. In this paper, I calculate semi- 
analytic expressions of the shapes of photon rings around black holes described by a new Kerr-like 
metric which is valid for all spins. I show that photon rings in this spacetime are affected by two 
types of deviations from the Kerr metric which can cause the ring shape to be highly asymmetric. 

I argue that the ring asymmetry is a direct measure of a potential violation of the no-hair theorem 
and that both types of deviations can be detected independently if the mass and distance of the black 
hole are known. In addition, I obtain approximate expressions of the diameters, displacements, and 
asymmetries of photon rings around Kerr and Kerr-like black holes. 

Subject headings: accretion, accretion disks - black hole physics - Galaxy: center - galaxies: individual 
(M87) - gravitation - gravitational leasing: strong 


1. INTRODUCTION 

The supermassive black holes at the centers of the 
Milky Way and M87 are the prime targets for high- 
resolution very-long baseline interferometric (VLBI) ob¬ 
servations with the Event Horizon Telescope (EHT), 
a planned global array of (sub)millimeter telescopes 
(Doeleman et al. 2009a,b; Fish et al. 2009). Recent 
VLBI observations at 230 GHz with a three-station ar¬ 
ray consisting of the James Clerk Maxwell Telescope 
(JCMT) and Sub-Millimeter Array (SMA) in Hawaii, the 
Submillimeter Telescope Observatory (SMTO) in Ari¬ 
zona, and several dishes of the Combined Array for Re¬ 
search in Millimeter-wave Astronomy (CARMA) in Cal¬ 
ifornia resolved structure on scales of only Ars (Sgr A*; 
Doeleman et al. 2008) and 5.5rs (M87; Doeleman et 
al. 2012), where rs = 2GM/c^ is the Schwarzschild ra¬ 
dius of a black hole with mass M and where G and c 
are the gravitational constant and the speed of light, re¬ 
spectively. Similar observations also detected time vari¬ 
ability on these scales in Sgr A* (Fish et al. 2011). Such 
measurements have demonstrated the feasibility of VLBI 
imaging of Sgr A* and the supermassive black hole in 
M87 on (sub-)event horizon scales. Simulations based 
on larger telescope arrays support the possibility of di¬ 
rectly imaging the shadows of these two supermassive 
black holes with the EHT (Doeleman et al. 2009a). 

Using the measurements of the mass and distance 
of Sgr A* from the monitoring of close stellar orbits 
(Ghez et al. 2008; Gillessen et al. 2009), Broderick et 
al. (2009a, 2011a) combined the VLBI observations of 
Sgr A* with measurements of its spectral energy distri¬ 
bution (Yuan et al. 2004 and references therein; Mar- 
rone 2006) and obtained values of the spin magnitude 


and direction employing a radiatively inefficient accre¬ 
tion flow (RIAF) model (Yuan et al. 2003; Broderick 
& Loeb 2006a). Likewise, Huang et al. (2009) ana¬ 
lyzed the VLBI data of Sgr A* using an accretion flow 
model with plasma wave heating for several different val¬ 
ues of the spin and disk inclination. Moscibrodzka et 
al. (2009) and Dexter et al. (2009, 2010) fitted the 
VLBI and spectral data to sets of images obtained from 
three-dimensional general-relativistic magnetohydrody¬ 
namic (3D-GRMHD) simulations (Gammie et al. 2003; 
Fragile et al. 2007, 2009; McKinney & Blandford 2009). 
In the future, the determination of the spin and orien¬ 
tation of the black hole can be further improved with a 
multiwavelength study of polarization (Broderick & Loeb 
2006b; see, also, Schnittman & Krolik 2009, 2010). 

The shadow is a prominent feature of resolved accre¬ 
tion flow images of supermassive black holes. The shape 
of the shadow depends uniquely on the mass, spin, and 
inclination of the black hole (e.g., Falcke et al. 2000) as 
well as on potential deviations from the Kerr metric (Jo¬ 
hannsen & Psaltis 2010b; hereafter JPlOb). Images of 
optically thin accretion flows also reveal a characteristic 
bright and narrow ring which surrounds the shadow, the 
so-called photon ring. This ring is the projection along 
null geodesics of photon trajectories that wind around 
the black hole many times. Thanks to the long path 
length through the emitting medium of the photons that 
comprise the ring, these photons can make a significantly 
larger contribution to the observed flux than individual 
photons outside of the ring. Images of the shadows of 
black holes with optically and geometrically thin accre¬ 
tion disks and of their photon rings have been calcu¬ 
lated by several authors (Bardeen 1973; Gunningham 
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& Bardeen 1973; Luminet 1979; Sikora 1979; Fukue & 
Yokoyama 1988; Jaroszynski et al. 1992; Viergutz 1993; 
Falcke et al. 2000; Takahashi 2004; Beckwith & Done 
2005; Bozza 2008; JPlOb; Chan et al. 2013). This pho¬ 
ton ring is clearly visible in several 3D-GRMHD simula¬ 
tions reported to date (Moscibrodzka et al. 2009; Dexter 
et al. 2009; Shcherbakov & Penna 2010). 

The photon rings of Sgr A* and the supermassive black 
hole in M87 have angular sizes of ^ 53 /ias and ^ 22 /ias, 
respectively, corresponding to a diameter of ^ 5.2rs (Jo- 
hannsen et al. 2012). Since the shape of the photon 
ring of a black hole is determined only by the geome¬ 
try of the underlying spacetime, it is independent of the 
complicated structure of the accretion flow making it an 
excellent target of VLBI imaging observations. For a 
Schwarzschild black hole, the photon ring is circular and 
centered on the black hole. For Kerr black holes with a 
spin 0 < |a| < GM/c^, where a = GJ/cM is the spin 
parameter of a black hole with angular momentum J, 
the photon ring retains a nearly circular shape and is 
displaced off center (Takahashi 2004; Beckwith & Done 
2005; JPlOb). The shape of the photon ring around Kerr 
black holes becomes asymmetric only in the case of ex¬ 
tremely high spin values a > 0.9GM/c^ if the inclination 
angle is large (JPlOb; Chan et al. 2013). Since the ring 
diameter of a Kerr black hole depends primarily on the 
black hole mass, imaging the ring around Sgr A* with 
VLBI can improve its current mass measurement in con¬ 
junction with the monitoring of the ephemerides of stars 
in its vicinity (Johannsen et al. 2012). 

Images of photon rings can be significantly altered 
if the black hole is not described by the Kerr metric 
(JPlOb). According to the no-hair theorem, the Kerr 
metric is the unique black hole metric in general relativ¬ 
ity and black holes are fully characterized by their masses 
and spins (see, e.g., Heusler 1996). However, while there 
exists evidence for the presence of event horizons in as- 
trophysical black holes (e.g., Narayan et al. 1997; Brod¬ 
erick et al. 2009b), the Kerr nature of black holes still 
remains untested to date. Several approaches have been 
developed to test the no-hair theorem which are based on 
Kerr-like metrics that deviate from the Kerr spacetime 
in parametric form (e.g., Manko & Novikov 1992; Collins 
& Hughes 2004; Clampedakis & Babak 2006; Vigeland 
& Hughes 2010; Johannsen & Psaltis 2011b; Vigeland 
et al. 2011). Observational signatures of potential viola¬ 
tions of the no-hair theorem can then be studied in either 
the electromagnetic (e.g., Johannsen & Psaltis 2010a,b, 
2011a, 2013; Bambi & Barausse 2011; Bambi 2012a,b, 
2013a,b; Krawczynski 2012) or gravitational wave spec¬ 
tra (see Cair et al. 2013; Yunes & Siemens 2013 for 
reviews). 

JPlOb analyzed the shapes of photon rings of compact 
objects described by a quasi-Kerr metric (Clampedakis & 
Babak 2006). This metric contains a quadrupole moment 
of the form 


Q = -M 



( 1 ) 


where the parameter e is dimensionless and measures po¬ 
tential deviations from the Kerr metric. When e = 0, 
the quasi-Kerr metric reduces to the familiar Kerr met¬ 
ric. For nonzero values of the parameter e, photon rings 


acquire an asymmetric shape. Therefore, the ring shape 
can be used as a direct measure of a potential violation 
of the no-hair theorem (JPlOb). 

By construction, however, the quasi-Kerr metric is only 
appropriate for the description of black holes that do not 
spin rapidly (Clampedakis & Babak 2006; Johannsen & 
Psaltis 2010a; Johannsen 2013a). In order to study de¬ 
viations from the Kerr metric for black holes with arbi¬ 
trary spins |a| < GM/c^, Johannsen & Psaltis (2011b) 
constructed a Kerr-like metric that can be used to model 
accretion flows even around rapidly-spinning black holes. 
Similar to the quasi-Kerr metric, this metric depends on 
one set of free parameters in addition to the mass and 
spin of the compact object. While this structure is ade¬ 
quate for the study of deviations from the Kerr metric of 
a certain type, it is not the most general form to param¬ 
eterize deviations from the Kerr metric, because such a 
Kerr-like metric should depend on at least four indepen¬ 
dent functions that measure these deviations instead of 
only one. Furthermore, as in the quasi-Kerr metric, the 
shapes of the shadows and photon rings of black holes in 
this metric can only be calculated numerically, because 
geodesic orbits are generally not integrable. Images of 
photon rings and accretion flows around black holes and 
some exotic objects in other theories of gravity were also 
analyzed by Bambi & Yoshida (2010), Amarilla et al. 
(2010), Bambi et al. (2012), Amarilla & Eiroa (2012), 
and Bambi (2013b). Chen & Jing (2012) and Liu et 
al. (2012a) studied the strong gravitational leasing near 
Kerr-like compact objects. 

In this paper, I study the shapes of photon rings around 
black holes that are described by a new Kerr-like metric 
proposed recently (Johannsen 2013b). This metric de¬ 
pends on four sets of deviation functions and describes 
a black hole for all values of the spin |a| < GM/c^. As 
in the case of the Kerr metric, the equations of motion 
of particles on geodesic orbits in this metric are fully 
integrable thanks to the existence of three independent 
constants of motion. This property allows for the cal¬ 
culation in closed form of the shapes of photon rings in 
the image plane of an observer at a large distance from 
the black hole in the limit of vanishing ring thickness. 
In this case, the photon ring reduces to the projection 
along null geodesics of the circular photon orbit in the 
equatorial plane of the black hole, which delineates the 
silhouette of the shadow of the black hole if its accretion 
flow is optically thin. 

I derive semi-analytic expressions of such photon rings 
around the Kerr-like black holes described by this met¬ 
ric and of the positions and shapes of these rings in the 
image plane. I show that the location of the circular 
photon orbit in the equatorial plane of the black hole as 
well as the shape of the photon ring depend on the mass, 
spin, and inclination of the black hole and on two types 
of deviations from the Kerr metric. In addition, I obtain 
approximate expressions of the ring diameter, displace¬ 
ment, and asymmetry in terms of these parameters for 
Kerr and Kerr-like black holes from least-squares fits of 
a large set of calculated photon ring images. 

In Section [2l I summarize the Kerr and Kerr-like met¬ 
rics that I use in this paper as well as some of their prop¬ 
erties. In Section[31 I determine the location of the circu¬ 
lar photon orbit and calculate expressions of the shapes 
of photon rings in the image plane. In Section|4l I derive 
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expressions of the diameters, displacements, and asym¬ 
metries of these photon rings and find approximate fit 
formulae. I discuss my conclusions in Section [5l From 
here on, I use geometric units and set G = c = 1. In these 
units, length (and time) scales are expressed in terms of 
the mass of the black hole, where M = GM/c^ = rs/2. 


2. KERR AND KERR-LIKE BLACK HOLE METRICS 
In Boyer-Lindquist coordinates, the Kerr metric is 


given by the elements 




2Mr\ 


K _ 
dtcj) 


S J 
2Mar sin^ 9 


g^=- 

arr ^ 5 




= s. 


= \r^ + 


2Mo?r sin^ 9 


sin^ 9, 


where 


A = — 2Mr -|- a^, 

E = r^ -I- cos^ 9. 


( 2 ) 


( 3 ) 


The Kerr-like black hole metric constructed by Jo- 
hannsen (2013b) is given by the elements 


gtt — ■ 


■ 


E[A — a^A 2 (r)^ sin^ 9] 

[(r^ -I- a?)Ai{r) — a^A^ir) sin^ 0]^ ’ 

a[(r^ -I- a^)Ai(r)A2(r) — A]S sin^ 9 
[{r'^ + aP‘)Ai{r) — a^A 2 {r)s\v? 9Y 

t 


9rr — 


AA5(r) 

3ee = E, 

E sin^ 9 [(r^ -|- a^)^Ai(r)^ — a^A sin^ d] 


500 = ■ 


[(r^ -I- a^)Ai(r) — a^A 2 (r) sin^ 0]^ 


where 




n—3 

CO 


A2{r) = l + Y,a2n(^) , 

n=2 


\ r 


^5(?') = l + ) , 

n=2 

E = E + /(r), 




( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


n—3 


The metric in Equation dl contains four free func¬ 
tions, /(r), Ai(r), A 2 (r), and A 5 (r), each of which de¬ 
pends on a set of parameters which measure potential 
deviations from the Kerr metric. The Kerr metric is re¬ 
covered when all deviation parameters vanish, i.e., when 


/(r) = 0, Ai(r) = A 2 (r) = A 5 (r) = 1. Thanks to 
these deviation functions, this metric can describe black 
holes in alternative theories of gravity that differ from 
Kerr black holes. However, instead of focusing on any 
particular modified gravity theory, which would a pri¬ 
ori limit the study of non-Kerr black holes to include 
only those that are predicted by that theory, it was de¬ 
signed in a phenomenological approach, which encom¬ 
passes large classes of modified theories of gravity corre¬ 
sponding to different choices of the deviation functions. 
Examples of such theories, which admit non-Kerr black 
holes that can be described by the metric in Equation (jH) 
in the appropriate limits, include Chern-Simons grav¬ 
ity and Einstein-Dilaton-Gauss-Bonnet gravity (see Jo- 
hannsen 2013b). Astrophysical observations can then be 
used to obtain constraints on deviations from the Kerr 
metric or to infer properties of the underlying, but so far 
unknown modified theory of gravity should the no-hair 
theorem be violated. 

Just as the Kerr metric, this metric is stationary, ax- 
isymmetric, and asymptotically flat. Since all station¬ 
ary, axisymmetric, asymptotically flat, vacuum metrics 
in general relativity can be written in terms of four in¬ 
dependent functions (see, e.g., Wald 1984 for a detailed 
discussion), it is natural to introduce four deviation func¬ 
tions to the Kerr metric instead of only one as in the 
quasi-Kerr metric (Glampedakis & Babak 2006) and the 
metric proposed by Johannsen & Psaltis (2011b). While 
the latter two metrics are thus less general than the met¬ 
ric in Equation (jd]) , they remain valid frameworks for the 
study of observational signatures of potential deviations 
from the Kerr metric. 

The deviation functions in Equations ([S])--®, in turn, 
are written as power series in M/r. Such a choice allows 
one to parameterize the deviations in a physically rea¬ 
sonable way given the enormous freedom of potential de¬ 
partures from the Kerr metric. The coefficients of these 
series are chosen so that the deviations they introduce 
are consistent with all current weak-field tests of general 
relativity such as those performed in the solar system 
and with double neutron stars (see Will 2006). Since one 
would expect that the lowest-order coefficients in these 
power series have the strongest effect on the observed 
signals, one could image to constrain these coefficients 
successively order by order with increasing accuracy of 
the measurements. 

For these reasons, I will only study black holes and 
their photon rings in this paper in the lowest-order met¬ 
ric, i.e, a given black hole only depends on its mass M, 
spin a, and the deviation parameters £ 3 , ai 3 , 022 , and 
052 . In the next section I will show that in this metric 
photon rings only depend on two deviation parameters, 
013 and 022 , which modify the (t,t), (t, <()), and 
components of the metric. Ultimately, constraints on the 
coefficients of higher order terms in the deviation func¬ 
tions will be required to obtain a deeper understanding 
of the full character of the deviation functions. 

While the deviation parameters £ 3 , 013 , 022 , and 052 
can be either positive or negative, they have to be larger 
than certain lower bounds, because otherwise the central 
object is a naked singularity and not a black hole. These 
bounds are given by the expressions (Johannsen 2013b) 

£3 > B3, 013 > 7?3, 
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where 


0:22 > B2, 052 > B2, 


B2 = - 

B3 = - 


(M +VM2-a2)^ 
(M +VM2-a2)^ 


( 10 ) 


( 11 ) 


Smaller values of these parameters are excluded. For 
all allowed values of the deviation parameters, the event 
horizon is independent of these parameters and coin¬ 
cides with the event horizon of a Kerr black hole, 


r+ = M + \/]VP - a2. (12) 

Thanks to the stationarity and axisymmetry of the 
metric in Equation (|4l), the energy E and angular mo¬ 
mentum Lz about the z-axis of particles moving along 
geodesics are conserved. The Kerr metric possesses a 
third conserved constant of motion, the so-called Carter 
constant (Carter 1968). Similarly, the metric in Equa¬ 
tion (j4|) admits a Carter constant of a similar form which 
is given by the expression (Johannsen 2013b) 

C = pg—{Lz—aE)'^+ii^a'^ cos^ 9-\ -^^(L^—aii'sin^ d)^, 

sin 0 

(13) 

where pg/T, = p^ is the 0 -component of the particle’s 
four-momentum and p is the rest mass of the particle. 

The existence of three constants of motion makes 
geodesic motion in both of these metrics fully integrable. 
This property greatly simplifies ray tracing algorithms in 
these spacetimes which are essential for the calculation 
of images of black hole accretion flows, because the equa¬ 
tions of motion of particles and, in particular, of photons 
can be written in first order form (see, e.g., Rauch & 
Blandford 1994; Dexter & Agol 2009 for such algorithms 
in the Kerr metric). 


3. CIRCULAR PHOTON ORBIT AND PHOTON RING 

In this section, 1 study the location of the circular 
photon orbit in the equatorial plane of a black hole de¬ 
scribed by the metric in Equation (U) and its projection 
along null geodesics onto the image plane of a distant 
observer. The image of the circular photon orbit marks 
the silhouette of the black hole shadow (Bardeen 1973; 
Cunningham & Bardeen 1973) which corresponds to the 
inner edge of the photon ring. Photon rings have a small 
thickness which is caused by the projection along null 
geodesics of photon paths that trace many orbits around 
the black hole but escape eventually. These curves are 
generally complicated functions of the radius r and the 
polar angle 9. In the equatorial plane, however, the tra¬ 
jectory is circular and located at a radius rph, the radius 
of the circular photon orbit. In this paper, 1 will ne¬ 
glect the thickness of the photon ring and only discuss 
thin photon rings which are the projection along null 
geodesics of the circular photon orbit. 

For a Schwarzschild black hole, the circular photon or¬ 
bit is located at the radius rph = 3M. In the Kerr metric, 
this radius is a function of the spin of the black hole. As 
the spin increases, the circular photon orbit approaches 
the event horizon and merges with the coordinate radius 


of the event horizon when a = M. In the Kerr-like met¬ 
ric given by Equation 0 , the location of the circular 
photon orbit also depends on the deviation parameters 
(c.f., Johannsen & Psaltis 2010a, 2011b). 

In order to calculate the location of the circular pho¬ 
ton orbit in the metric given by Equation 0 , I make 
use of the expressions of the energy E and axial angular 
momentum Lz of massive particles on circular equatorial 
orbits which I already calculated in Johannsen (2013b). 
In the limit of vanishing rest mass, the energy and axial 
angular momentum diverge at the location of the cir¬ 
cular photon orbit. Both the energy and axial angular 
momentum for massive particles on circular equatorial 
orbits (including the ISCO, the innermost stable circular 
orbit) generally depend on the spin a and the deviation 
parameters ea, aia, and a 22 but are independent of the 
deviation parameter q; 52 . The latter deviation parame¬ 
ter only affects the (r, r) component of the metric, which 
plays no role for particles on circular equatorial orbits 
in this regard (Johannsen 2013b). The location of the 
circular photon orbit, however, depends only on the spin 
and the deviation parameters ais and q; 22 , because all 
terms that are affected by the parameter 63 also contain 
the rest mass of the particle which vanishes for photons. 
I, then, calculate the radius of the circular photon orbit 
by numerically finding the root of the inverse of the en¬ 
ergy E. Equivalently, the axial angular momentum Lz 
could be used instead, which would yield the same result. 

In Figure [TJ I plot the radius of the circular photon 
orbit as a function of the spin of the black hole for differ¬ 
ent values of the parameters ai 3 and q; 22 - In Figure [H 
I plot contours of constant circular photon orbit radius 
as a function of the spin and the deviation parameters 
ai 3 and q; 22 - The dependence of the circular photon 
orbit on these deviation parameters is similar to the de¬ 
pendence of the ISCO on the same deviation parameters 
(see Johannsen 2013b). The location of the circular pho¬ 
ton orbit depends strongly on both deviation parameters 
ai 3 and a 22 except for black holes that have a nearly 
maximal spin a ^ M. In that case, the circular pho¬ 
ton orbit radius depends only weakly on the deviation 
parameters. The location of the circular photon orbit 
decreases for increasing values of the spin and decreasing 
values of the parameter q;i 3 . The circular photon orbit 
radius also decreases for increasing values of the param¬ 
eter a 22 if a > 0 and decreasing values of the parameter 
a 22 if a < 0. For nonspinning black holes, the radius of 
the circular photon orbit is independent of the parameter 
Oi22- 

In the following, I derive semi-analytic expressions of 
the images of photon rings around black holes that are 
described by the metric in Equation (jd]). In this deriva¬ 
tion, I make use of the existence of three independent 
constants of motion in this metric. 

First I dehne an inertial reference frame of an observer 
at a large distance from the black hole (or of any local 
observer). For the basis vectors, I choose 


et = Cet +760, 
1 



1 


^9 ! - 
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-1 -0.5 0 0.5 1 
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Fig. 1. — Radius of the circular photon orbit as a function of 
the spin for different values of the deviation parameters (top) ais 
and (bottom) 0:22 • In each panel, only one deviation parameter is 
varied, while the other one is set to zero. Top panel: the radius 
of the circular photon orbit decreases for increasing values of the 
spin and decreasing values of the deviation parameter 0 : 13 . For 
negative values of the parameter 013 , the circular photon orbit 
merges with the event horizon when a = M. Bottom panel: the 
radius of the circular photon orbit decreases for decreasing values of 
the parameter 022 if a < 0 and increasing values of the parameter 
0^22 if a > 0. If a = 0, the circular photon orbit is located at the 
radius r = 3M, irrespective of the value of the parameter 022 - 



where {et, e,., eg, e^} is the coordinate basis of the metric. 
The basis vectors of the local observer are related to the 
coordinate basis of the metric via the relations 

e& = e'bc;,, (15) 



-1 -0.5 0 0.5 1 

Spin (M) 



-1 -0.5 0 0.5 1 

Spin (M) 


Fig. 2. — Contours of constant circular photon orbit radius as 
a function of spin and the deviation parameters (top) 0:13 and 
(bottom) 022 - In each panel, only one deviation parameter is varied 
and the other one is set to zero. In the top panel, the radius of the 
circular photon orbit decreases for increasing values of the spin and 
decreasing values of the parameter 0 : 13 . In the bottom panel, the 
radius of the circular photon orbit decreases for increasing values 
of the spin and decreasing values of the parameter 0:22 if a < 0 
and increasing values of the parameter 0:22 if a > 0. If a = 0, the 
radius of the circular photon orbit is independent of the parameter 
a 22 and is located at a radius of 3M. The dashed line corresponds 
to a Kerr black hole. The shaded region marks the excluded part 
of the parameter space. 

where 

= V&0 ( 16 ) 

and where = diag(—1,1,1,1) is the Minkowski met¬ 
ric. The constants ( and 7 in Equation (HI have to be 
chosen in a manner so that the local basis is orthonor¬ 
mal. From the relations in Equation (jl6l) . 1 obtain the 
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expressions 


c=, 


7 = - 


9<t><t> 




g<t><t> 


gti, 


g<p<p y 5(0 gttg<i><f> 


(17) 


in terms of the metric elements given by Equation O- 
In the new basis, the locally measured energy and axial 
angular momentum are given by the expressions 


P* = CE - 


1 


p'^ = 


^/g<t><t> 


-L,. 


(18) 

(19) 


I can then define the impact parameters x' and y' in the 
usual manner (Bardeen 1973), which are the Cartesian 
coordinate axes of the image plane of the observer located 
at a distance r = tq and inclination 0 = i: 


X =-ro^ = -ro 

pt 


/_ / 

y =rQ^ = ro- 

pt 


(i + gj«) ’ 

±\/ 0 (O 


v»c (> + IS«) 

where (see Johannsen 2013b) 


/ = 


Pe _ ±\/Q(») 


yfg^ 

Q{i) = r] + cos^ i — cot^ i. 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


In these expressions, all metric elements are evaluated at 
r = tq and 9 = i and I introduced the parameters 


^ E ’ 
C 


v= 


E 2 ’ 


(24) 

(25) 


where C is the generalized Carter constant given in Equa¬ 
tion (fT^ . 

In the limit tq —>■ oo, I obtain the expressions 


X = -- 




Sin 2 


y' = ±^J<cl{i), 


(26) 

(27) 


which are similar to the ones for the Kerr metric 
(Bardeen 1973). Note that the y'-axis is chosen so that 
it lies in a plane with the rotation axis of the black hole 
and that the origin of the image plane is centered on the 
black hole. 

In order to obtain explicit expressions of the param¬ 
eters ^ and rj which correspond to the photon ring of 
the black hole, I determine these parameters from the 
location of the circular photon orbit. Since these param¬ 
eters are combinations of the constants of motion, they 
are conserved along the null geodesics that project the 
circular photon orbit onto the photon ring in the image 
plane and can be calculated at any point along these 


geodesics. This is most easily achieved at the location of 
the circular photon orbit. 

At the circular photon orbit, the r-component of the 
photon four-momentum (Johannsen 2013b), 

f = kAtiM, (28) 

s 

where 


i?(r) = [{r^W)A,{r)-aA 2 {r)e-m-a)"+v]: (29) 


as well as its radial derivative vanish. Since from Equa¬ 
tions ©-(III]) S > 0 and A 5 (r) > 0, I only need to 
consider the function i?(r). Thus, I simultaneously solve 
the equations 


i?(r) = 0, 
dR{r) 


(30) 


for the parameters ^ and rj. The expressions for these 
parameters are given explicitly in the Appendix. 



x’ (M) 


Fig. 3.— Images of photon rings around a Kerr black hole with 
spin a = 0.998M for different values of the inclination i. Photon 
rings are shown at inclination angles 2 = 0° (dark blue), i = 10° 
(navy blue), i = 90° (red). Even for a Kerr black hole with 
such an extreme spin, the ring shape is nearly circular except for 
high inclinations. 


Images of photon rings around Kerr black holes have 
already been calculated by several authors (e.g., Bardeen 
1973; Luminet 1979; Takahashi 2004; Beckwith & Done 
2005; JPlOb; Chan et al. 2013). As an example, I plot 
in Figure [3] images of photon rings around a Kerr black 
hole with spin a = 0.998M for different values of the 
observer inclination. As discussed in Takahashi (2004), 
JPlOb, and Chan et al. (2013), these images are nearly 
circular except for very large values of the spin a > 0.9M 
and high inclination angles. The image size is primarily 
determined by the mass of the black hole. For increasing 
values of the spin at a fixed inclination f > 0 °, the image 
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Fig. 4.— Photon rings around black holes at an inclination angle i = 90° for different values of the spin a and deviation parameters 0:13 
(top row) and 0.22 (bottom row). The ring size is largely determined by the black hole mass and the parameter 0 : 13 , while it is affected only 
marginally by the spin and the parameter OL 22 • The displacement is primarily affected by the spin of the black hole and depends weakly on 
the deviation parameters. Nonzero values of the deviation parameters, however, can cause the photon rings to be significantly asymmetric. 


is displaced off center of the image plane. If the black 
hole is non-rotating (a = 0 ) or viewed face-on (i = 0 °), 
the photon ring is centered on the black hole. The dis¬ 
placement is reminiscent of the location of the caustics 
in the Kerr spacetime (Rauch & Blandford 1994; Bozza 
2008). 

In Figure m I plot images of the photon rings of black 
holes that are described by the Kerr-like metric in Equa¬ 
tion (|4|) for different values of the spin and the deviation 
parameters ai 3 and a 22 at an inclination angle i = 90°. 
The size of the photon ring is altered primarily by the 
parameter ai 3 and increases for increasing values of this 
parameter, while the parameter 022 affects the image size 
only marginally. Both deviation parameters only have a 
weak impact on the image displacement compared to the 
effect of the spin unless the spin is very large, but they 
have a strong impact on the shape of the image which 
can be significantly distorted and asymmetric. However, 
for positive values of the parameter 013 as well as for neg¬ 
ative values of the parameter a 22 , photon rings around 
Kerr-like black holes can be less asymmetric than photon 
rings around Kerr black holes with the same spin if the 
spin is very large. Nonetheless, significant ring asymme¬ 
tries can only be caused by nonzero deviations from the 
Kerr metric unless the spin is near-extremal (a > 0.9M) 
and the inclination is large. 


4. DIAMETER, DISPLACEMENT, ASYMMETRY 

Following Takahashi (2004) and JPlOb, I define the 
(horizontal) displacement D of the photon ring relative 
to the center of the image plane by the expression 


p _ l^max ^minl 

2 

(31) 

where and a;' ■ are the maximum and minimum 

abscissae of the ring, respectively. Thanks to the reflec¬ 
tion symmetry across the equatorial plane of the black 
hole and the choice of the coordinate system the image 
plane, there is no displacement in the vertical direction. 
Following JPlOb, I further define the average radius of 
the ring by the expression 

1 

Ric. 

(32) 

where 


R=yJ{x' - Df 

(33) 

y' 

tana= —. 
x' 

(34) 

The ring diameter is then 


L = 2{R). 

(35) 
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Fig. 5. — Diameter (left), displacement (center) and asymmetry (right) of the photon ring around a Kerr black hole as a function of the 
inclination for different values of the spin a in the range a = 0 (dark blue), a = O.IM (navy blue),..., a = 0.9M (orange), a = 0.998M 
(red). The diameter is practically constant at a given value of the spin and depends only very weakly on the observer inclination. The 
displacement is determined by both the spin and the inclination angle. The asymmetry, however, is small except for very large values of 
the spin a > 0.9M at large inclination angles. 


The ring asymmetry is defined by the expression 


{R-{R))"da 

V 27r 


(36) 


Takahashi (2004) calculated the ring displacement as 
a function of the spin for several inclination angles and 
determined approximately the maximum ring diameter 
(across the geometric center of the ring) from a set of 
calculated ring images. JPlOb calculated the ring di¬ 
ameter, displacement, and asymmetry for a number of 
images generated by a ray tracing algorithm in the Kerr 
and quasi-Kerr spacetimes. They obtained approximate 
expressions for these ring characteristics in terms of the 
mass, spin, and quadrupolar deviation parameter e which 
was defined in Equation O- Chan et al. (2013) re¬ 
fined the expressions for the diameter and asymmetry of 
the photon ring around a Kerr black hole using a set of 
photon ring images simulated by their high-performance 
ray tracing code GRay. In the following, I derive semi- 
analytic expressions of the ring diameter, displacement, 
and asymmetry using the expressions of the ring image 
that I obtained in the previous section. 

Since the ring coordinates x' and y' in the image de¬ 
fined in Equations (1^ and (EZl) are parameterized in 
terms of the radius r, it is convenient to use Equa¬ 
tion (1341) in order to convert the integrals in the ring 
diameter and asymmetry over the angle a into integrals 
over the radius by the coordinate transformation 

da = H{r)dr, (37) 


where 


H{r) 


d arctan 

x' 


fdy' tdx' 

dr y dr 

dr 


xl2 _|_ yl2 


(38) 


is the Jacobian of this transformation. The ring diameter 
and asymmetry then become: 

2 

L=— H{r)dr, (39) 

Jri 


A = 2\ 


{R-{R)) H{r)dr 


(40) 


where ri and r 2 are the two roots of the discriminant in 
the y'-coordinate in Equation (|27l) with ri <r 2 - Note the 
extra factor of two in the integration over radius, since 
the integrals have to be evaluated twice in this radial 
interval. 

In Figure [SJ I plot the diameter, displacement, and 
asymmetry of the photon ring around a Kerr black hole 
as a function of the inclination for various values of the 
spin. Chan et al. (2013) obtained very similar results 
from fits to a number of simulated photon ring images. 
For small values of the spin, the diameter is practically 
constant with a value ~ 10.4M (see JPlOb). For values 
of the spin 0 < |a| < 0.77M, the ring diameter lies in 
the range lOM ^ L < 10.4M. For large values of the 
spin a > 0.77M, the ring diameter can be smaller than 
lOM. The displacement depends approximately linearly 
on the spin and the sine of the inclination angle (JPlOb; 
Chan et al. 2013; c.f., Bozza 2008). The asymmetry is 
generally small and only significant at very large values 
of the spin a > 0.9M and high inclinations. 

In Figure [H I plot the diameter, displacement, and 
asymmetry of the photon ring around a black hole de¬ 
scribed by the metric in Equation (|4]) with a spin a = 
0.9M as a function of the inclination for different values 
of the deviation parameters 0:13 and 022 - The ring di¬ 
ameter depends only weakly on the parameter 0:22 and is 
practically constant for fixed values of the parameter ai 3 
ranging from L Ri 8M for ai 3 = —2 to L Ri 10. 8 M for 
ai 3 = 2. Therefore, the value of the deviation param¬ 
eter Q!i 3 could be inferred from a measurement of the 
ring diameter if the mass and distance to the black hole 
are known with sufficient accuracy. The ring displace¬ 
ment, on the other hand, is determined primarily by the 
spin of the black hole. Even at an inclination angle of 
i = 90°, values of the deviation parameters 0:13 = ±2 or 
a 22 = ±2 change the displacement caused by the spin 
alone by at most 50%. While the spin could, in princi¬ 
ple, be inferred from the ring displacement, it is unlikely 
to yield a useful spin measurement in practice, because 
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Fig. 6.— Diameter (left), displacement (center), and asymmetry (right) of photon rings around Kerr-like black holes with a spin a = 0.9M 
as a function of the inclination angle for different values of the deviation parameters ai3 (top row) and 0:22 (bottom row). The ring diameter 
depends only weakly on the parameter 022, while it is practically constant for fixed values of the parameter 013. The deviation parameter 
013 cs,n be inferred from a measurement of the ring diameter if the mass and distance to the black hole are known. The displacement of 
the ring is altered by the parameters 013 and 022 but depends primarily on the spin. Negative values of the parameter 013 and positive 
values of the parameter 022 can cause the ring shape to be significantly more asymmetric than the photon ring around a Kerr black hole 
with the same spin. The ring asymmetry is a direct measure of the degree to which the no-hair theorem is violated. 


the exact location of the black hole in the image would 
have to be determined. The asymmetry of the ring is 
very large if ais <C 0 or if a 22 ^ 0. For values of the de¬ 
viation parameters ai 3 = —2 and 022 = 2, respectively, 
the asymmetry is ~ 2.7 and ~ 3 times as large as the 
asymmetry of the ring around a Kerr black hole with the 
same spin at an inclination oi i = 90°. As in the case of 
the quasi-Kerr metric (see JPlOb), the ring asymmetry 
is, therefore, a direct measure of the degree to which the 
no-hair theorem is violated. 

The modified displacement and asymmetry of pho¬ 
ton rings around Kerr-like black holes correspond to the 
shifted location of the circular photon orbit as shown in 
Figure [TJ For negative values of the parameter ai 3 as 
well as for positive values of the parameter a 22 , the ra¬ 
dius of the circular photon orbit decreases if the black 
hole spin lies in the range 0 < a < M. Therefore, the 
photons that comprise the image of the circular photon 
orbit will experience stronger lightbending which causes 
the image to become more displaced and more asym¬ 
metric as in the case of a the ring around a Kerr black 
hole with the same spin. Likewise, at least for deviations 
that are measured in terms of the parameters ai 3 , the 
ring diameter decreases. 

Next, I obtain approximate expressions of the diame¬ 


ters, displacements, and asymmetries from a large num¬ 
ber of simulated photon rings. First, I obtain fit for¬ 
mulae for these properties of photon ring around Kerr 
black holes. I computed the diameters, displacements, 
and asymmetries of 17,113 photon rings from Equations 
dMl), dSB), and go]) choosing spin values in the range 
from a = O.OIM to a = 0.99M in steps of Aa = O.OIM 
and from a = 0.991M to a = 0.999M in steps of 
Aa = O.OOIM for values of the inclination ranging from 
i = 0.01 to i = 1.57 in steps of Ai = 0.01. 

I make the following ansatz for the diameter, displace¬ 
ment, and asymmetry of the photon ring, where the sub- 
and superscripts “K” refer to a Kerr black hole: 


LK(ba)- 

s Lf -1- Lf cos {nfi + , 

(41) 

DKii,a)K 

s Df sin (nfi) , 

(42) 

AK(i,a)s 

1 

(43) 

“ Afi^f -b Afi< ' 


Using this ansatz, I perform a least-squares fit and obtain 
the expressions 

Lf = 10.3907 - 0.520134a2-^®®^'^ - 0.109423a^^-®'^°®, 
Lf = 0.141644a2 “®^'‘ -f 0.077576a^-®°®^® 
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A 


K 

1 


n 


K 

3 


A 


K 

2 


nf 


-0.305764a'‘‘‘^®®^ + 0.603316a®-®®^^'‘ 
-0.801615a®®™^ + 0.892182a^®-^®®® 
-0.508699a^®•®l°^ 

-2.00026 - 0.0424525a®-^^“^ 
+0.0503142a^®-®^^^ + 0.118102a^^^-^®°, 
3.14246 + 0.0970934a®'‘°®^2 
+0.0964881a^^-^®^^ - 0.0990121a^®® ®®®, 

2 a + 0.188324a^-2®°2® + 0.153158a^°-^'‘^® 
+0.0982701a'‘®-^®25 0.05008470^2® °®^, 

1 - 0.0149883a° °®°^°2 + 0.0188647a°-^29558 
+0.0370128a2'4°°®° + 0.147758a®°-2°25 
-0.163190a®^-^°^‘‘ - 0.0148721a2®°-^'‘2^ 

3.29826 4.06676 1.69824 


(1 _ q)0.0708700 ' q1.96447 ' ^2.06984 

+0.535821a®®°°®° + 0.801151a^°®®^®, 
-2.55018 + 0.679338(1 - a)° ®°°®®‘‘ 

+0.2632290^-^2755 o.ll3085o2 ®^°®4 

-0.0657765o^°-°°'‘® - 0.0497249o°° ®°®2, 
0.0482563 1.17042 


(1 - o) 


0.291648 


,1.99594 


- 0.2596720^-®°®®^ 


-0.306290o®°-^°°2 _ o.323404a^^®-®2°, 

0.235657 + 1.39539(1 - a)°-^®®^26 

+1.334680^-^2230 _ 0.8424820^-®°°®® 

+0.09627170®°-°®®® + 0.140634o®®°-®°® 
+0.1272510°°°-"^®®. (44) 


In the fit of the displacement, 1 held the first terms 
in the parameters Di and n 2 fixed so that D^{i,a) ~ 
2osini as found in JPlOb. The fit formula of the ring di¬ 
ameter is accurate to < 0.09% for spin values a < 0.997M 
and to < 0.16% for spins o < 0.999M. The fit formula 
of the displacement is accurate to < 1 . 1 % for spin val¬ 
ues a < 0.97M, to < 2.4% for spins a < 0.98M, and 
to < 3.7% for spins a < 0.999M. For values of the 
asymmetry A > O.OIM, the asymmetry fit is accurate 
to < 1% for spin values o < 0.6M, to < 4.8% for spins 
a < 0.98M, and to < 19.4% for spins a < 0.999M. The 
largest uncertainties in the fit occur at lo-w inclinations, 
and the fit is accurate at all spin values to < 2.5% for 
inclinations i > 22.4°. Values of the ring asymmetry 
A < O.OIM, -which occur only at spin values a ~ 0 and 
inclinations i ~ 0 , were partly affected by numerical un¬ 
certainty and, therefore, neglected in the fit. Since these 
values of the asymmetry are very small, this affects the fit 
only marginally. The asymmetry can also be fitted with 
an ansatz of the form Aq sin" i as in JPlOb and Chan et 
al. (2013). Such a fit, however, introduces comparatively 
large errors at high spins, where the asymmetry deviates 
significantly from a sinusoidal form. 

Next, 1 obtain approximate expressions for the diam¬ 
eters, displacements, and asymmetries of photon rings 
around the Kerr-like black holes described by Equa¬ 
tion ®. 1 calculated the diameters, displacements, 

and asymmetries of 288,300 photon rings with spins 
a/M = 0.05, 0.10,..., 0.95, 0.998, inclinations i = 
0.1, 0.2,..., 1.5, and deviations ai 3 =—1, —0.9,..., 2, 


022 = —1, — 0.9,..., 2. For the respective fit functions, 

1 introduced polynomial functions of the deviations pa¬ 
rameters 013 and 022 to the corresponding fits for Kerr 
photon rings in Equations (I41l) -"(l43 )l given by the expres¬ 
sions 

L{i, a, oi3, 022 ) ~ + T 1 O 13 + T 2 O 13 + L^a 22 

+Lia%2 + ^^(1 + 7 ^ 5013 ) 
x(l + L 6 a 22 ) cos (nfi + , (45) 

D{i, a, oi3, 022 ) ~ of' (1 + Diais + D 2 a^ 3 ) 

X (1 + 743 O 22 + 7 I 4 O 22 ) 

X sin [nf{l + niai 3 + 772013 ) 

X (1 + 773022 + 774022 ) 7 ] , (46) 

A(i, a, 013 , 022 ) ~ (1 + A 1 O 13 + A 2 O 13 + ^ 3 *^ 13 ) 

X (1 + A 4 O 22 + A^(X22 4” 4 I 3 O 22 ) 

X [A]^(1 + A70i3)(1 + ^§ 022 )*"^ 

+^^(1 + ^9<ai3)(l + ^ 10022 )*"^] 

(47) 

From these, 1 obtained the following values for the fit 
parameters: 

Li = 0.390575 + 0.289090a®-°®®®° + 0.3514850^°-°°°^ 

L 2 = -0.0210469 - 0.0704973a®-°°°®® - 0.156894a®°-°°®°, 
L 3 = -0.0771662a®-®°°°° - 0.112522a®°-®®®®, 

Li = -0.0232177a°-®'‘^'‘® - 0.0474855a'‘°®-°'‘®, 

Ls = -0.0910695 + 0.5978610®"^-°®®°, 

L 6 = -0.665419al-'^^'‘°^ 

Di = -0.0471772 - 0.0714876a'‘-°’'®^® - 2.80254o®°^-°°\ 
D 2 = 0.0100617 + 0.1973160^^-°®^^ + 0.1736940®°“^-®°°, 

D 3 = 0.175643 + 0.107683a°-°'‘°°° - 0.309073a®°-'‘®®'‘, 

Di = 0.0144201 + 0.3231820®°-®°®^ - 0.377931a^®^-®°°, 

771 = -0.0617871 - 0.151930a®^-°®®^ + 2.21857a®®®-°°°, 

772 = 0.00753226 - 0.00323579a®^-®°®^ + 1.40633a°^-°°°°, 

773 = 0.0180080 - 0.0970060a2®-°‘^®i + 0.4398630^°^-^®°, 

774 = -0.0183142 - 0.2281960°®-^^^° + 0.0583285a°°-°®^®, 
Al = 0.0403176 - 0.458487a®-^'^°^'‘, 

A 2 = -0.0649604 + 0.0847185a®-®°°®°, 

A 3 = 0.00650092 - 0.00620249o“°-°°^ 

A 4 = 1.24363 + 0.211407a°°-™^°, 

As = 0.281333 - 0.00508372a°-°°°°^ 

Ag = 0.0287109 - 0.266727a®-®°®®°, 

At = 0.286452 + 0.640543a°^-°®^®, 

As = 0.014822 - 0.407110a°-°°'^'‘®, 

Ag = 0.274555 - 0.722911a®-°®®°°, 

Aio = 0.208396 + 0.799466a°°-^^°®. (48) 

The fit of the diameter is accurate to < 7.5% for spins 
a < 0.95M and to < 28.5% for spins a < 0.998M. The 
fit of the displacement is accurate to < 13.3% in spin 
range O.IM < a < 0.85M with an average accuracy of 
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2%. The accuracy is significantly smaller at high and 
very low spins and the error can exceed 100 % in some 
cases. Finally, the fit of the asymmetry has an average 
accuracy of < 4.6%, but it can deviate from the simu¬ 
lated data set by factors of order unity, especially at very 
high spins. If 013 = 022 = 0, these fits reduce to the cor¬ 
responding fits for photon rings around Kerr black holes 
in Equations (|4T]) - (I^ . 

5. CONCLUSIONS 

In this paper, I analyzed the location of the circular 
photon orbit and the shapes of photon rings around Kerr- 
like black holes which are described by a recently pro¬ 
posed metric (Johannsen 2013b). This metric depends on 
four independent parameters and can serve as a frame¬ 
work to model accretion flows around black holes with 
arbitrary spin |o| < M, for which observational signa¬ 
tures of potential violations of the no-hair theorem can 
then be studied (c.f., Johannsen & Psaltis 2010a). This 
work extends the previous analysis of photon rings in 
JPIOb who studied potential violations of the no-hair 
theorem in terms of a quasi-Kerr metric (Glampedakis 
& Babak 2006) up to values of the spin |a| < 0.4M. I 
showed that the radius of the circular photon orbit de¬ 
pends on only two deviation parameters which can shift 
its location significantly. I derived semi-analytic expres¬ 
sions of the photon ring in the image plane as observed 
at a large distance away from the black hole and of the 
inferred diameter, displacement, and asymmetry of the 
ring. From large sets of simulated ring images, I cal¬ 
culated approximate formulae of the ring diameter, dis¬ 
placement, and asymmetry of both Kerr and Kerr-like 
black holes. I determined the accuracy of these formulae 
which is very high for Kerr black holes and moderate for 
Kerr-like black holes. 

In the case of a Kerr black hole, the ring diameter de¬ 
pends only very weakly on the spin and the inclination 
and can, therefore, be used to improve the mass measure¬ 
ment of Sgr A* in combination with the monitoring of 
stars on close orbits around the supermassive black hole 
(Johannsen et al. 2012). For photon rings in the new 
Kerr-like metric, I showed that the image size depends 
primarily on the mass and the deviation parameter 013 . 
Therefore, this parameter could be determined from the 
image size if the mass and distance are known with suf¬ 
ficient accuracy. Currently, Sgr A* has a measured mass 
and distance with uncertainties of only Ri 9% and r: 5%, 
respectively (Gillessen et al. 2009), which are compara¬ 
ble to the changes of the ring diameter introduced by the 
parameter 013 . The uncertainty of the mass of the su¬ 
permassive black hole at the center of M87 is currently 
reported to be R 29% assuming no uncertainty in the 
distance (see Giiltekin et al. 2009 and references therein; 
see, however, Walsh et al. 2013). 

As pointed out by JPIOb, the displacement and asym¬ 
metry of the photon ring only need to be measured in 
units of the ring diameter in order to infer the black 
hole spin and deviation from the Kerr metric. An ab¬ 
solute measurement of these quantities, which would be 
affected by the uncertainties in the mass and distance 
measurements, is not necessary. While the displacement 
of the ring depends primarily on the spin of the black 
hole, a measurement of the spin from the displacement 
will be very difficult in practice, because the exact posi¬ 


tion of the black hole in the image would have to be de¬ 
termined. The asymmetry, on the other hand, depends 
strongly on the deviation parameters. Since photon rings 
around Kerr black holes are nearly circular except for 
values of the spin a > 0.9M and at large inclinations 
(JPIOb; Chan et al. 2013), the ring asymmetry is a di¬ 
rect probe of a violation of the no-hair theorem (c.f., 
JPIOb). 

The EHT is expected to achieve an unprecedented an¬ 
gular resolution which lies significantly below the angular 
sizes of the photon rings of Sgr A* and the supermas¬ 
sive black hole in M87. Thanks to the planned inclusion 
of the Atacama Large Millimeter/submillimeter Array 
(ALMA) and state-of-the-art receivers at the other sta¬ 
tions, the EHT will also be highly sensitive. Especially 
the measurement of closure phases between the stations 
will allow for a highly-resolved study of the shadows and 
the photon rings of these sources (c.f., Broderick et al. 
2011b). Johannsen et al. (2012) estimated the flux den¬ 
sity of the photon ring of Sgr A* to be on the order of 
0.1 Jy which means that the ring should be sufficiently 
bright to be imaged by the EHT. 

In order to assess the prospects of the EHT to detect 
potential deviations from the Kerr metric, detailed sim¬ 
ulations of the photon ring will be required which take 
into account the actual source flux as well as the reso¬ 
lution and sensitivity of the telescope array. A recent 
analysis of the current VLBl data in the context of a 
RIAF accretion flow model in a quasi-Kerr background 
placed first constraints on such deviations (Broderick et 
al. 2013). Depending on the magnitude of the deviation 
parameter, these models also predict different amounts 
of correlated flux density which should be distinguish¬ 
able at least along baselines between ALMA and sites in 
North America (see Figure 5 in Fish et al. 2013). 

Here 1 limited my analysis to include only thin photon 
rings which delineate the shadow silhouette of a black 
hole if its accretion flow is optically thin, which is the 
case at least for Sgr A* at submillimeter wavelengths. 
Actual photon rings, however, have a small thickness and 
the ring images 1 calculated in this paper coincide with 
the inner edges of these photon rings, which are marked 
by a sharp flux fall-off. Further refined models of the ring 
which take into account its thickness will be required for 
an image analysis, but the tools developed here should be 
very useful for such models. Due to the relativistic effects 
of boosting and beaming, the ring thickness is not uni¬ 
form and is largest near the equatorial plane on the side 
of the black hole that approaches the observer. Kam- 
ruddin & Dexter (2013) recently constructed a simple 
crescent model to incorporate the varying ring thickness 
(see, also, L. Benkevitch et ah, in preparation). 

In practice, the image of the photon ring is slightly 
blurred by the array resolution and interstellar scattering 
and needs to be disentangled from the image of the accre¬ 
tion flow in the ww-plane, which is only partially sampled 
by a given set of baselines. Johannsen et al. (2012) used 
the location of the nulls of the visibility function as a 
benchmark for the uncertainties in such a measurement 
and argued that it should be feasible with the EHT. A 
full analysis will need to extract the ring with a pattern 
matching technique of the full visibility function. 

In the case of Sgr A*, the no-hair theorem may also be 
tested with two different techniques in the regime of weak 
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gravitational fields. The presence of the spin of the black 
hole as well as its quadrupole moment introduces torques 
on stars in close orbits around Sgr A* which cause their 
orbital planes to precess (Will 2008). While this effect 
is too small for the currently known S-stars in the vicin¬ 
ity of Sgr A*, because their average distances from the 
center of mass are too large, these precession frequencies 
might be measurable for stars within ^ lOOOrg of the 
black hole, if they are monitored over a sufficiently long 
time (Merritt et al. 2010). For these stars, potentially 
perturbing effects such as those of other stars (Sadeghian 
& Will 2011) and dark matter (Sadeghian et al. 2013), 
as well as drag forces exerted by the hot plasma sur¬ 
rounding Sgr A* (Psaltis 2012), may be negligible. Such 
observations could be obtained using the second genera¬ 
tion instrument GRAVITY for the Very Large Telescope 
Interferometer (Eisenhauer et al. 2011). 

In addition, if one of the stars orbiting around Sgr A* is 
a pulsar, timing observations may measure characteristic 


spin-orbit residuals that are induced by the quadrupole 
moment of the black hole (Wex & Kopeikin 1999; Pfahl 
& Loeb 2004; Liu et al. 2012b). Recently, Mori et al. 
(2013) reported the discovery of the first pulsar in the 
vicinity of the Galactic center which lies at a distance 
of only ^ 0.1 pc (Rea et al. 2013) corresponding 
to ^ 2.5 X lO^rg. Future measurements obtained by 
the monitoring of stars, timing observations of pulsars, 
and VLBl imaging are largely independent of each other 
and will be affected by different systematic uncertainties. 
The combination of their results could lead to a robust 
test of the no-hair theorem. 
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APPENDIX 


The set of Equations (l30l) has two solutions. One of them corresponds to one single point in the image plane, which 
is unphysical. Therefore, this solution has to be rejected. The parameters ^ and i] in Equations are then 

given by the expressions 


c=- 


—3a‘*M^ai3 -I- a^r(7M‘*ai3 — 5M^rai3 + Mr^ + r'^) — 2M^rai3 — + r'^) 

ar[M2a22(2a^ -I- 3r2) — 5M^ra22 — Mr^ + 

r{a®A'/’^ai3Q;22(Yf^ai3Q;22 + 6Mr^Q!i3 


^ o2r®[M2a22(2o^ -I- 3r^) — 5M^ra22 — Mr^ + r^]'^ 

—4r^a22) — M^r[2M'^ a'l^a22 + ^^fa\^a22 + ^‘2-M^r'^a\^a22 — 2M^r^ ai^a22{ctiz + 4 ^ 22 ) 

-|-8M^r"‘Q;i3(a22 “ 3ai3) -I- M^r®ai3(15ai3 -I- 34a22) + 6Mr®a22(ai3 - 2 q; 22) - 8r’^a22] 
-|-a'*M^r^{4M®a^3Q;22 — M^ra\^a\2 + SM'^r^al^a22 + M®r^ai3 [8022 — ai3(6a22 + 49)] 
-|-2M®r"‘Q;i3[38Q;i3 -|- Q;22(2a22 + 35)] - + 8Q;i3a22 + 5 ^ 22 ) + 2M^r®( 222 ( 0:13 - O 22 ) 

(23(222 ~ 22(213) + 2Mr®(5ai3 + 8022 ) + 8 r®a 22 } + — 2M^^al^al2 + 

-|-4M®r^(2i3(222 - 2M’^r^(2i3 [ 013(022 + 21) + 8 ( 222 ] + 2M®r"‘(2i3 [ 28(213 + 022(4022 + 15)] 
-M^r^ai3{19ai3 + 30a22) + 2MV® [ 013 ( 4(222 + 21) - 160 ^ 2 ] " [24ai3 + (15 - 8022 ) 022 ] 

-|-2M^r®(7oi3 — 2022 ) + 10 Mr®a 22 + 4r^°} — r^°(3AT^oi3 — 2 M^roi 3 — 3Mr^ -P r^)^}. 

In the Kerr limit, these parameters reduce to the familiar expressions (Bardeen 1973) 

r^(r — 3M) -P a^(r -P M) 


? = - 


1 = 


a{r — M) 
r^[4a^M — r(r — 3M)^] 
a?{r — MY 


( 1 ) 


( 2 ) 

( 3 ) 

( 4 ) 
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